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Abstract 

The finite-temperature magnetic properties of Fe x Pdi_ x and Co x Pti_ x alloys have been inves- 
tigated. It is shown that the temperature-dependent magnetic behaviour of alloys, composed of 
originally magnetic and non-magnetic elements, cannot be described properly unless the coupling 
between magnetic moments at magnetic atoms (Fe,Co) mediated through the interactions with in- 
duced magnetic moments of non-magnetic atoms (Pd,Pt) is included. A scheme for the calculation 
of the Curie temperature (lb) for this type of systems is presented which is based on the extended 
Heisenberg Hamiltonian with the appropriate exchange parameters Jy obtained from ab-initio 
electronic structure calculations. Within the present study the KKR Green's function method has 
been used to calculate the Jij parameters. A comparison of the obtained Curie temperatures for 
Fe x Pdi_a; and Co x Pt\- x alloys with experimental data shows rather good agreement. 

PACS numbers: 71.20.Be, 75.30.Hx, 75.40. Cx 
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I. INTRODUCTION 



Whether a magnetic material is technologically useful or not depends on its properties 
at finite temperatures. However, the ab-initio treatment of finite-temperature magnetism 
remains a challenge despite ongoing progress in this field during the last decades. 

In this context, itinerant-electron 3d-transition metals and their alloys receive particu- 
lar interest. For these systems, finite-temperature magnetic properties cannot be described 
successfully neither within the collective-electron Stoner model nor within the local-moment 
model based on the Heisenberg Hamiltonian (for an overview see Refs. l|, |2|). The Stoner 
model, which treats the transition to the paramagnetic state as vanishing of the local mag- 
netic moments, accounting thereby only for the longitudinal spin fluctuations, grossly overes- 
timates the Curie temperature Tc- A bigger success has been achieved using the Heisenberg 
model, which accounts for temperature-induced transverse spin fluctuations and character- 
izes the paramagnetic state by orientational disorder within the system of localized magnetic 
moments. However, a magnitude of the moments is assumed to be unchanged upon fluc- 
tuations. In some recent studies the Heisenberg model approach has been combined with 



ab-initio band-structure ca 
eters from first principles 



culations, that allow to evaluate the exchange coupling param- 
3)]. In this way, trends of the Curie or Neel temperature with 
composition can be quantitatively described for many systems, in particular, for transition 
metal monoxides |4|, dilute magnetic semiconductors |5| or transition metals [6|. 

Despite the rather satisfying results obtained within this combined approach for several 
systems, the description of finite-temperature magnetism of transition metals and alloys still 
suffers from many problems owing to the restrictions of the Heisenberg model. For some 
itinerant-electron systems, e.g. Ni, the thermally-induced longitudinal spin fluctuations 
play a crucial role in describing properly the temperature-dependent magnetisation and 
obtaining the correct value for the critical temperature. A phenomenological theory of 
finite temperature magnetism, which accounts for both types of fluctuations on the same 
footing, was developed in the past jl- 1Q|. This theory was used in combination with ab-initio 
electronic structure calculations to describe temperature dependent magnetic properties of 
Fe, Co and Ni llMl3|. In this way a much better agreement with experiment, as compared to 



calculations based on the Heisenberg model, was obtained. In particular, proper accounting 
for longitudinal fluctuations results in the vanishing of local magnetic moments on Ni atoms 
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above Tq 13J, in agreement with experiment. 

Other interesting itinerant-electron systems in this context are alloys or compounds com- 
posed of originally magnetic and non-magnetic elements. Such systems exhibit the so called 



covalent magnetism 
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15] , where magnetisation of the 'non-magnetic' atoms is governed 



by the spontaneously magnetised atoms via the strong spin- dependent hybridisation of their 
electronic states. The Fe x Pdi_ x and Co x Pti_ x alloys considered in the present work belong 
to this type of systems. To describe the temperature-dependent magnetism of such systems 
on the basis of Heisenberg model, one obviously has to account properly for the behavior 
of the Pd/Pt sub-lattices. Only a few ab-initio studies of finite-temperature magnetism of 
systems of this kind have been done so far. Similarly to the work mentioned above jllL 



13J these studies were based on a generalisation of the classical Heisenberg Hamiltonian in 



one or another way to account for the different character of magnetism on different types 
of atoms. Mryasov et al. [16| have investigated the compound FePt and showed that the 
anomalous temperature-dependence of its magneto-crystalline anisotropy energy (MAE) is 
due to the induced Pt magnetic moments. In another study of these authors 17| , a crucial 
role of the magnetic moment induced on Rh was demonstrated for the stabilisation of the 
ferromagnetic state of FeRh and for the control of the antiferromagnet-ferromagnet phase 



transition. Lezaic et al. 



18| emphasized the need to account for longitudinal fluctuations 



of magnetic moments induced on Ni atoms for a proper description of the temperature- 
dependence of the spin-polarisation at the Fermi energy Ep in the half-metallic ferromagnet 
NiMnSb. Sandratskii et al. 19] investigated several ways for accounting for the induced 
magnetic moments within the spin-spiral approach used for the calculation of the exchange 
coupling parameters in NiMnSb and MnAs. Using these results Sandratskii et al. 19) have 
studied finite-temperature magnetic properties of NiMnSb; their findings are consistent with 
the findings of Lezaic et al [18| . 

In this work we introduce an ab-initio method to describe finite-temperature magnetism 
of systems with spontaneous and induced magnetic moments. The method is based on an 
extension of the Heisenberg Hamiltonian by adding a term which describes the induced 
magnetic moments within the linear response formalism. Our approach relies on a combi- 
nation of ab-initio band-structure calculations with Monte Carlo (MC) simulations based 
on the extended Heisenberg model. To test this approach, we investigate finite-temperature 
magnetic properties of Pd-rich Fe^Pdi.^ alloys, with Fe concentrations up to 20 at.%, as 
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well as of ordered and disordered CosPt, CoPt and CoPt3 alloys that are interesting both 
for fundamental reasons and for possible use in industrial applications because of their high 
magnetic anisotropy 



20 



21] . Theoretical investigations of finite-temperature magnetism of 



these systems failed so far to reproduce experimental results with satisfactory accuracy [22 ] . 
We demonstrate in the present work that a combination of ab-initio band-structure calcula- 
tions with the Monte Carlo (MC) simulations based on the extended Heisenberg model gives 
satisfying agreement between theoretical and experimental values of critical temperatures. 
We found that despite their small magnitudes, the moments induced on non-magnetic atoms 
(Pd, Pt) have an important influence on finite-temperature magnetic order. 



II. EXPERIMENTAL DETAILES 



The Fe^Pdi-a; films were thermally evaporated onto oxidized silicon substrates from sepa- 
rate effusion cells for Pd and Fe in an ultra high vacuum system (base pressure 510 -11 mbar). 
The film thicknesses were between 15 and 20 nm. The deposition rate of the two components 
could be controlled independently, resulting in an accuracy of the Fe concentration of 1%. 
Auger spectroscopy on a sample with nominally 7% Fe content provided an independent 
value of 7.8% Fe for this film. After deposition, the films were patterned into a six-terminal 
Hall-bar geometry. Measurements of the anomalous Hall effect at temperatures between 2K 
and 300 K provided the magnetization M(T), from which the Curie temperatures of the 
films were deduced. 



III. THEORETICAL APPROACH 
A. Ground-state calculations 



Within the present work, spin-polarised electronic structure calculations for the ground- 
state have been performed using the multiple scattering KKR (Korringa-Kohn-Rostoker) 



Green's function method [23j] in the scalar-relativistic approximation. The local spin density 
approximation (LSDA) for density functional theory was used with the parametrisation for 



the exchange-correlation potential due to Vosko, Wilk, and Nusair ! 



24j. The potential was 



treated within the atomic sphere approximation (ASA) with the radii of the spheres around 
Fe/Co and Pd/Pt sites chosen by requiring the ratios of the corresponding volumes to be the 



same as for the pure elements. For the angular momentum expansion of the Green's function 
a cutoff of l max = 3 was applied. For substitutional^ disordered alloys, the self-consistent 
coherent potential approximation (CPA) method was employed. A geometry optimisation 
was performed, i.e., the lattice constants of the alloys have been obtained by minimisation 
of the total energy. 

B. Extended Heisenberg Hamiltonian 

The finite temperature properties of the investigated systems were studied by Monte 
Carlo simulations based on the Heisenberg model with the underlying Hamiltonian given by 

He* = "E Jt M MA - £ 

ij ij 

~ E J"™*™! ■ (!) 

ij 

Here the classical Hamiltonian was generalised to allow an application to itinerant- 
electron systems consisting of magnetic and non-magnetic atoms having magnetic moments 
Mi and rrij, respectively, connected with corresponding exchange coupling parameters. The 
dependence of the induced magnetic moments rrii on a specific magnetic configuration is 
treated via linear response formalism (for details see the Appendix). 

We suppose that the induced magnetic moments on Pd or Pt atoms are governed only by 
the magnetic moments of the Fe or Co atoms in Fe^Pdi-^ and Co-jPti-.,; alloys, respectively, 
arranged in the first neighbor shell around the non-magnetic atom, so that 

= E xr M Mj = xr M E 4- • (2) 

jeM jeM 
The notation J^ gM means that the sum includes only such terms where j is a site with an 
inducing magnetic moment. 

In this case (see Appendix) the susceptibility X^~ M can be approximated by the value 
found for the ground state of a system with well defined collinear spontaneous (Fe, Co) and 
induced (Pd, Pt) magnetic moments: 



The exchange coupling parameters Jja between atoms % and j in the Eq. ([I]) were obtained 
via the formula of Lichtenstein et al. [3] : 

i r^F 

J tJ = -— Im J dE Tr - t^) 7* - f-i) r f . (4) 
with the relation 



Jij 



J i 



\ m i\\ m j\ 

In Eq. PJ, t icr and are the spin (a) and site dependent single-site and scattering 
path operator matrices occurring within the KKR formalism 25]. 



C. Evaluation of Tq 

I — I 

The Curie temperature Tc was evaluated with the Monte Carlo (MC) method [26[ us- 
ing the standard Metropolis importance sampling algorithm [271 ]. on the basis of the model 
Hamiltonian in Eq. (pQ). The number of atoms in the MC unit cell for different concen- 
trations was taken between 1728 and 4000 and Tq was determined from the peak position 
of the temperature-dependent susceptibility. The Fe/Co magnetic moments were treated 
during MC simulation as localised and changed only their orientation. On the other hand, 
the magnetic moments on Pd/Pt atoms could change their absolute value as well as the 
orientation in accordance with the changed magnetic configuration around these atoms. Eq. 
([2]) implies that the magnetic moments on Pd/Pt are proportional to the vector sum of mag- 
netic moments at neighboring Fe/Co atoms with only nearest neighbors taken into account. 
This means that each MC step consists of 1) change of the orientation of a magnetic moment 
on the Fe/Co atoms, 2) search for all nearest neighbour Pd/Pt atoms and calculation of the 
orientation and absolute value of their moments using the susceptibilities X m ~ M via Eq. (|3]). 
The change of the energy of the total system is due to both effects. For disordered alloys 
the resulting Tq values in addition were averaged over up to 20 different configurations. 

Obviously, the approach described above, accounts for the contribution of spin polarised 
'non-magnetic' atoms to the exchange interactions between the magnetic atoms. As this 
contribution is temperature-dependent, it allows a corresponding description of the temper- 
ature dependent magnetisation. In particular, it accounts for longitudinal spin fluctuations 
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occurring on the 'non-magnetic' sub-lattice. Below it will be shown that this rather simple 
scheme gives rather good agreement with experimental data for the systems under consid- 
eration. 

IV. RESULTS FOR FE^PDi^ ALLOYS 
A. Ab-initio calculations 

The scheme for calculation of temperature dependent magnetic properties, described 
above, was used to investigate disordered Fe x Pdi_ z alloys with Fe concentration up to 20 
at.%. The exchange coupling parameters J Fe ~ Fe and J? e ~ Pd , shown in Fig. (TJ have a 
similar dependency on the distance Rij for all investgated alloys. In the Pd reach limit 
(Fe concentration x < 0.2) the exchange coupling parameters J Fe ~ Fe corresponding to the 
average distance between magnetic atoms is rather small and does not allow to create long- 
range magnetic order in the system, as was demonstrated by corresponding restricted MC 
simulations. 




FIG. 1: Exchange coupling parameters J Fe ~ Fe (a) and J Fe ~ Pd (b) for Fe x Pdi_ x alloys at different 
concentrations. The dashed lines represent results for ordered FePd3 (note that in this case for the 
Fe atoms the nearest neighbour interaction J Fe ~ Fe is absent because there are only Pd nearest 
neighbours) . 
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FIG. 2: Results of spin-spiral calculations for ordered FePd3 and disordered Feo.25Pdo.75: (a) the 
energy of spin spirals as a function of the wave vector q = ~(0,0,q z ), obtained for the disordered 
alloy (full circles), for the ordered compound with non-zero Pd magnetic moments (full squares) 
and for the ordered compound with Pd magnetic moments equal to (open squares); (b) - magnetic 
moment of inequivalent Pd atoms in ordered FePd3 as a function of the wave vector q z . Fe and 
Pd2 occupy the sites (0,0,0) and (^, |,0), respectively, Pdl atoms occupy the sites (0, ^, ^) and 

(i,o i). 



According to the experimental findings 



28], the alloys exhibit ferromagnetic order for 



e concentrations above 0.1%. In addition, previous experimental 



29 



and theoretical 



15 



37 



41] investigations on the magnetic properties of Fe^Pdi-^. alloys have shown a strong 
host polarisation by magnetic Fe impurities leading to a giant magnetic moment per impurity 
atom, up to 12.9^5. A wide-spread regime of magnetised Pd atoms leads to ferromagnetic 
order in diluted Fe^Pdi.^ alloys despite a large distance between the magnetic Fe atoms. 

The crucial role of the induced magnetic moment on Pd for the Fe-Fe exchange interac- 
tions can be demonstrated by an analysis of the energy of spin spirals as a function of a wave 
vector, shown in Fig. [2] for the disordered alloy Feo.25Pdo.75 m comparison with the results 
for the ordered compound FePd3. The calculations have been performed for spin spirals 
along the z direction with the Fe magnetic moments tilted by 90° with respect to the z-axis. 
For the ordered system an increase of the wave vector of the spin spirals is accompanied first 
by an increase in energy reflecting the stability of the ferromagnetic order in the system. A 
further increase of the wave vector above q z > n /2a, leads to a decrease of the energy of the 
spin spiral (Fig. [2^). This behavior is governed by a decrease of the Pd magnetic moments 



S 



at these wave vectors (see Fig. [2b), that diminishes their role in the Fe-Fe exchange. 

The role of Pd becomes clearly visible for the spin spirals in FePd3 with the Fe-Pd 
exchange interactions being suppressed. This can be achieved by forcing the Pd induced 
magnetic moments to be perpendicular to the Fe magnetic moments, and therefore to be 
equal to (see Appendix). In this case the minimum of the spin-spiral energy corresponds 
to an AFM state, i.e. at q = ir/a, that originates from Fe-Fe exchange interaction (see Fig. 
[Q). In the case of the disordered alloy, the dependence of the spin-spiral energy on the wave 
vector is different because the random distribution of the Fe atoms allows Fe atoms to be 
nearest neighbours with a strong FM interaction. Due to this, the system retains the FM 
order at all values of wave vector q. 
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FIG. 3: a) Pd DOS in Fe x Pdi_ x for x = 0.003 and x = 0.15 for spin-up (upper panel) and spin- 
down (lower panel) states. Ground state characteristics of Fe x Pdi- x alloys vs. Fe concentration: 
b) density of states of Pd at the Fermi level; c) Pd spin magnetic moments; d) Fe spin magnetic 
moments. 

While the ground-state magnetic properties of Fe^Pdi.^ alloys can essentially be under- 
stood on the basis of ab-initio electronic structure calculations, a theoretical description 
of finite-temperature properties faces many difficulties due to the itinerant-electron nature 
of magnetism. Theoretical investigations based on the results of ab-initio electronic struc- 
ture calculations have been performed for examp 
the model approach formulated by Bloch et al. 



e by Mohn and Schwarz [15|. They used 
to describe the magnetic behavior of 



a system characterised by the coexistent local- and itinerant-electron magnetism. Within 
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this approach the system is characterised by two interacting subsystems: (i) one having 
local magnetic moments showing a Curie- Weiss-like behaviour and (ii) an itinerant elec- 
tron subsystem magnetically polarised by the effective Weiss field, with the corresponding 
parameters found by ab-initio electronic structure calculations. 

The spin moment on every Pd atom is induced by the magnetic moment of the Fe atoms 
and all surrounding induced Pd magnetic moments (see Appendix). The rather big absolute 
value and large region of Pd induced magnetic moments around an Fe atom is a result of the 
high magnetic spin susceptibility of pure Pd and Fe^Pdi-^. alloys with small Fe concentration, 
that is determined by a large Pd density of states (DOS) at the Fermi level, n(Ep) (see Fig. 
Ek). m turn, n(Ep) decreases with the increase of Fe content in the alloy (Fig. [3b), resulting 
in a decrease of the partial magnetic susceptibility of the Pd atoms. Thus, at very small Fe 
concentrations the induced Pd spin moment can extend to big distances - inducing shell- 
by-shell a spin moment in the Pd subsystem. This polarisation mechanism decays with the 
distance from the magnetic impurity. When the Fe concentration increases, the regions with 



the induced moments overlap and, as was pointed in M. Shimizu et al. 43|, |44j, when the 
Fe concentration is larger than 0.1% the Pd magnetic properties can be described well by 
band calculations, as done in our present calculations using the CPA alloy theory. 

Figure [3l c and d, represent the spin magnetic moments of Pd and Fe versus the Fe 
content in Fe^Pdi-j;. As is seen, the Fe magnetic moments change only slowly with the 
increase of Fe concentration, while the variation of the Pd spin magnetic moments is rather 
pronounced. This can be directly connected to the decrease of the Pd DOS at the Fermi 
level. 



For a more detailed analysis, we investigated the properties of the induced Pd magnetic 
moment using the ab-initio calculations. In particular, we studied the distribution of the Pd 
magnetic moment in a Pd host in the limit of very small Fe concentrations, i.e. around a 
single Fe impurity. To see that the induced magnetic moment at every Pd atom is determined 
not only by the Fe magnetic moment but also by the surrounding Pd magnetic moments, one 
can compare the un-enhanced induced Pd spin moments created by only one Fe atom with 
the total induced spin moments in Pd. The total induced magnetic moment distribution in 
Pd can be found by solving the system of equations (1A3|) within a selected region around an 
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R 0j /a 

FIG. 4: Magnetic moment distribution in Pd around a single Fe impurity, as a function of the 
distance from the magnetic atom. Full squares represent the Pd magnetic moments self-consistently 
obtained within a cluster with 3 atomic shells of Pd around an Fe impurity embedded into a Pd 
host, full circles give the results obtained for the same system but with the exchange potential 
switched off. The open circles represent the induced magnetic moments in Pd calculated within 
linear response formalism using Eq. (jAljl in the Appendix. 

Fe atom (see Appendix). In addition, in the present work the moment distribution has been 
obtained by self-consistent electronic structure calculations instead of using linear response 
formalism. Fig. H] shows the slow decay of the induced magnetic moment with the distance 
(full squares). The corresponding un-enhanced induced moment in these calculations has 
been obtained by suppressing the effective exchange B-field for the Pd atoms during the 
SCF-cycle. These un-enhanced magnetic moments compare very well with those obtained 
from linear response formalism (Eq. (1A1H ) (open circles). These are shown in Fig. H]also 
for larger distances. One can see that the decrease with the distance of the Pd un-enhanced 
spin moment is very fast compared to the enhanced one. For the nearest Pd neighbors of 
Fe atom these values differ approximately by a factor of 2, while the difference for the next 
nearest neighbors is already an order of magnitude. The local exchange enhancement, well 
approximated within the linear approach at small values of the induced magnetic moments, 
should keep the ratio of these two values approximately constant. The obtained results give 
evidence for a more complicated picture of the creation of the induced magnetic moment in 
accordance to the description given in the Appendix. 

The effect of temperature-induced magnetic disorder within the Fe-subsystem was anal- 
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ysed within ab-initio calculations, describing magnetic disorder within the uncompensated 
Disordered Local Moment (DLM) approximation. Using this approximation an effective 
alloy of two types of Fe atoms with opposite spin directions and having different concen- 
trations is treated using the CPA alloy theory. In this way one can study the dependence 
of the induced magnetic moment of individual Pd atoms on the average magnetic moment 
in the system. Fig. [5] shows that the magnetic disorder in the Fe subsystem (assumed to 
be temperature induced) is accompanied by a decease of the total magnetic moment in the 
system and results in a decease of the induced magnetic moment in the Pd subsystem. One 
can see a rather good linear dependence of the induced Pd magnetic moment as a function 
of the magnetic moment of the Fe subsystem, for nearly all Fe concentrations. Only in the 
limiting case of low Fe concentration (1%), a noteworthy deviation from linear behavior is 
observed. This deviation will influence the final results in a Curie temperature evaluation 
correspondingly. 




FIG. 5: Results of ab-initio calculations for induced Pd moment in three Fe-Pd alloys, using the 
uncompensated DLM approximation as a function of the Fe average magnetic moment. 



B. Finite-temperature magnetism of Fe^Pdi x alloys 

The temperature dependent magnetic properties of Fe^Pdi-a, alloys were investigated by 
performing Monte Carlo simulations. They show an absence of an ordered FM state for 
the alloys with Fe concentration up to x ~ 0.2 if the Fe-Fe exchange interactions mediated 
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by Fe-Pd interactions are neglected. This clearly demonstrates the importance of these 
interactions. The Curie temperature for disordered Feo.2Pdo.8 alloy in this case is around 
60K, much lower than observed experimentally (around 400 K). To illustrate the role of 
Fe-Pd interactions in the formation of magnetic order, Fig. [H] shows spin configurations 
obtained within MC simulations for (T = 0.1K) for the Fe x Pdi_ x alloy with x = 0.02. 
Fig. represents a result for the case that the Fe-Pd exchange interactions are neglected 
and therefore it shows only the magnetic moments of Fe. Fig. [6b represents the spin 
configuration when the Fe-Pd first-neighbour interactions are taken into account and shows 
only those Fe and Pd magnetic moments which give a contribution to the total energy of a 
system given by Eq. ([T]). 

A comparison of Tq obtained within the MC simulation based on the Hamiltonian in Eq. 
([I]) with the experimental data for FexPdi-^ alloys is shown in Fig. [7J Obviously, a rather 
good agreement is obtained for the whole concentration range. It should be emphasized 
once more that all parameters for the model Hamiltonian (Eq. (p])) are obtained within 
ab-initio electronic structure calculations using the scheme described in the Appendix. Of 
course, going beyond the various approximations the final results can be improved to get 
better agreement with the experimental results. Fig. [JJ shows that the theoretical result 
obtained by Mohn and Schwarz 15| for Fe a Pdi_a; alloys at small Fe concentrations is also 
in a good agreement with experiment. However, it should be emphasized that this work is 
based on a semi-phenomenological approach. 




FIG. 6: The distribution of spin moments in unit cell used within the MC simulations without 
(a) and with (b) the induced Pd moments taken into account. Large arrows correspond to the Fe 
magnetic moments, small arrows to Pd magnetic moments. 
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FIG. 7: (a) The Curie temperature for different Fe concentrations in Fe x Pdi_ x alloys: present 
results vs theoretical results of Mohn and Schwarz [lj| and experimental data: Expt. 1 - present 



work, Expt. 2 
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system 


a 


^spin (Co) 


/^spin (Pt) 




a.u. 


fJ-B 




Co 3 Pt 


ordered 6.98 


1.82 


0.36 




disordered 7.06 


1.88 


0.25 


CoPt 


ordered 7.23 


1.93 


0.37 




disordered 7.23 


2.03 


0.27 


CoPt 3 


ordered 7.31 


1.76 


0.26 




disordered 7.37 


2.19 


0.25 



TABLE I: Equilibrium lattice constants and magnetic moments at Co and Pt atoms for ordered 
and disordered Co^Pti-a; alloys. 

V. RESULTS FOR CO.PTi^ ALLOYS 

A. Magnetic moments and exchange coupling constants 

The calculated equilibrium lattice constants for ordered and disordered Co^Pti-^ alloys 
are shown in Tab. [H For ordered CoPt, we used a simplified Llo geometry, assuming c=a 
instead of c=0.98a found in experiment. The magnetic moments of the Co and Pt atoms 
for the equilibrium lattice constants obtained in the scalar-relativistic mode are presented 
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in Tab 

study 



em More details about ground-state properties of Co-Pt can be found in our earlier 

The exchange coupling constants for the investigated systems were evaluated via Eq. 
01]). The dependence of on the distance between the atoms i and j is displayed in Fig. 
IH] with the left panels showing the situation when both i and j are Co atoms and the right 
panels showing the situation when i is a Co atom and j is a Pt atom. If experimental lattice 
constants were used instead of equilibrium lattice constants, the Jy constants would change 
slightly but both the trends and the values would remain similar as in Fig. [HJ 

One can see from Fig. [S] that for CoaPt and CoPt the coupling between the moments 
on the Co atoms do not differ very much from the results for their disordered counter 
parts, Coo.75Pto.25 and Coo.50Pto.50 respectively. However, the situation changes dramatically 
for CoPt3 (lower left panel in Fig. [8]). The pronounced difference between the data for 
the ordered and the disordered system stems mainly from the fact that there are no Co 
atoms present for some coordination spheres around a central Co atom in ordered CoPt3. 
Concerning the coupling between moments on Co and Pt atoms, the degree of long-range 
order has a larger influence than for the Co-Co coupling. For ordered alloys, the Jy constants 
significantly vary also with composition. For disordered alloys, on the other hand, the Jy 
constants do not vary very much with composition. For ordered CoPt3, there is a surprisingly 
strong Co-Co coupling between atoms which are 2.83 a apart. We verified that for larger 
distances no comparable strong coupling occurs. 



B. Curie temperatures 

The Curie temperatures Tq of the investigated Co-Pt systems evaluated by means of the 
Monte-Carlo technique are shown in Tab. [Til In addition, results for fee Co are given in this 
table. The two theoretical Curie temperatures correspond to two different Hamiltonians 
used to describe the magnetic coupling. The first T c (denoted as "Co-Co only" in Tab. 
ITT|) corresponds to the standard Heisenberg Hamiltonian for magnetic moments only on the 
Co atoms. The second Tc (denoted as "Co-Co and Co-Pt") corresponds to the extended 
Heisenberg Hamiltonian (Eq. (Jij), accounting for the coupling between moments on Co 
atoms Jfj°~ Co as well as for the coupling between moments on Co and Pt atoms Jy°~ pt , 
with the moments on Pt atoms determined via Eqs. (j2J) and ([3]). The values of Tc obtained 
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system model 


T c [K] T c [K] 




ordered disordered 



fee Co 


MC 




1100 






experiment 




1388 






Mohn-Wohlfarth 




3523 




Co 3 Pt 


MC, Co-Co only 




800 


750 




MC, Co-Co and Co- 


-Pt 


900 


880 




experiment 






1100 




Mohn-Wohlfarth 




1803 


1120 


CoPt 


MC, Co-Co only 




360 


620 




MC, Co-Co and Co- 


-Pt 


620 


760 




experiment 




727 


830 




Mohn-Wohlfarth 




1964 


850 


CoPt 3 


MC, Co-Co only 




-180 


370 




MC, Co-Co and Co- 


-Pt 


150 


520 




experiment 




288 


468 




Mohn-Wohlfarth 




241 


510 



TABLE II: Curie temperatures Tc for fee Co and for ordered and disordered CoxP\,\- x alloys. 



Experimental results 



47|, 



481 ] are shown together with results of our Monte-Carlo calculations with 
either both Co-Co and Co-Pt coupling or with only Co-Co coupling included. For comparison, 
results obtained by relying on the Mohn-Wohlfarth theory are also shown 
Tc implies antiferromagnetic ordering. 
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49 



50J] . Negative 



earlier 
et al. 



relying on the semi-empirical Mohn-Wohlfarth theory 5l| were taken from Qi 



49[ for fee Co, from Kashyap et al. [22| for ordered Co-Pt compounds and from 



Ghosh et al. [50| for disordered Co^Pti-a; alloys. Experimental values |47j, |48| are shown 
for comparison. Note that an experimental value for Tc for ordered CosPt is not available 
because the ordered phase is not stable for this composition. 

n 

As is seen, the Mohn-Wohlfarth theory [51J gives reasonable agreement with experiment 
at small Co concentration. However, increasing the Co content leads to large discrepancies 
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between the theory and experiment. This results from the limitations of the Mohn-Wohlfarth 
theory: it was developed for homogeneous itinerant-electron systems, which is not the case 
for Co^Ptx^ alloys. 

Our ab-initio scheme based on an extended Heisenberg Hamiltonian (both Co-Co and 
Co-Pt coupling included) accounts quantitatively for the trends of Tq with the composition 
and with the degree of long-range order. If the coupling mediated via moments at Pt atoms 
is not included, the results are unrealistic. This is especially true for ordered CoPt 3 , where 
an antiferromagnetic order is established at finite temperatures (reflected by a negative value 
for T c ) if the coupling between moments on Co and Pt is ignored. 

VI. CONCLUSIONS 

As was shown in the present work by the examples of Fe x Pdi_ x and Co^Pti-^ alloys, the 
finite-temperature magnetism of alloys composed of magnetic and non-magnetic elements 
requires to account for the exchange interactions between magnetic atoms, mediated by 
the exchange interaction with non-magnetic atoms. This inplies in particular that one has 
to account properly for the induced magnetic moment within the Monte Carlo simulations 
which are based in the present work on a corresponding extension of the standard Heisenberg 
Hamiltonian. The approach presented suggests to describe the induced magnetic moment 
on non-magnetic atoms within linear response formalism, being proportional to the vector 
sum of magnetic moments of neighboring magnetic atoms. This ansatz allows for substantial 
technical simplifications and leads to substantial improvement of the results when compared 
to simpler schemes. 

The finite-temperature calculations for Fe x Pdi_ x and Co :r Pti„ :E alloys performed within 
this approach give the dependence of Tq on the composition as well as on the degree of 
long-range order in good agreement with experimental data. The case of ordered CoPt 3 also 
demonstrates that even if the coupling between nearest inducing moments is antiferromag- 
netic, the magnetic order can still be ferromagnetic due to the effect of coupling between 
inducing and induced moments. A mere inspection of the Jij constants thus cannot serve 
as a reliable indicator of ferromagnetic or antiferromagnetic order at T 7^ OK. 
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Appendix A: Induced magnetic moments 

The magnetic moment induced on site i (Pd or Pt) by the exchange field Bj° due to 
magnetic moments at site j can be calculated within the linear response formalism using 



the expression 



52|: 



mf(f) = Im j dETra 

7T J 

x I d 3 r'G(r, f ', E)H* c (r')G(r\ r, E) 
= [ X ^r,f')B] c (r')d 3 r' (Al) 

Here Hj C (f') = aBj C (f') = aeMBj C (f') with a - the matrix of Pauli matrices [5^, e*M 
the unit vector in the direction of spontaneous magnetic moment of atom j, Bj c (f) the local 
exchange field at the site j. Note that, neglecting relativistic effects, the magnetic moment 

induced by a neighboring magnetic atom is parallel to the direction cm of the magnetic 
moment Mj of this atom. 

The total induced magnetic moments on a Pd or Pt atom is represented as a response to 
the exchange field of all surrounding atoms by the following expression: 



rrii{r) = 



E / XJT^ ?')Z7 M ^')d*r' (A2) 
+Xi [rjBi (r) . 



Analogous to J2jeM defined at Eq. (|2J), the sum J2j em means summation over sites with 
induced magnetic moments. 
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This can be reformulated in terms of local magnetic moments Mj of Fe (Co) and rrii of 
Pd (Pt) 



irrii 



(A3) 



with 



Mi= d 3 rMi(f); rhi = / d 3 rrhi(r) . (A4) 

Jflws J^ws 

Here we use the following reformulation for the first term in Eq. (1A2j) . that is more 

convenient for the model implementation: 



u M 



X ™- M (r,r')Bp M (r')d 3 r' (A5) 



TyXC,M / l\ 



For the second term, using a linearised expression for the exchange potential in the case 



of a small induced magnetic moment on Pd and Pt sites 
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54j , one can write analogously: 



/ X ™- m (r,r>)Bp m (r>)d 3 ri (A6) 



Solving the system of equations ( 1A3I) for a restricted region around a magnetic impurity 
atom gives the distribution of the induced magnetic moments on the non-magnetic atoms. 
This can be done for the ground state (T = OK). Alternatively, without any approximations, 
one can get these values within ab-initio calculations for embedded magnetic atoms by using 
the CPA alloy theory, assuming an uniform distribution of the induced magnetic moment. 
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Strictly spoken, one can go beyond the linear approximation in the expansion of the 
exchange potential. However, the linear approximation makes the use of this scheme in 
subsequent Monte Carlo simulations much easier. 

By making an additional simplification one can restrict to one response function X™~ M 
within the Monte-Carlo simulations. This quantity is defined to give the induced magnetic 
moment as a response to the exchange fields of only the surrounding nearest neighbour 
magnetic atoms: 
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FIG. 8: Exchange coupling constants between moments of a Co atom on site i and on other 
surrounding Co atoms (j) (left panels) and between moments on a Co atom (i) and on surrounding 
Pt atoms (j) (right panel), for ordered and disordered Co^Pti-z alloys. The horizontal axis shows 
the distance Rij between i and j atoms in units of lattice constants. 
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